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We observed the nucleus of M33, the third-brightest galaxy in the Local Group, 
with the Space Telescope Imaging Spectrograph at a resolution at least a factor 
O ' of 10 higher than previously obtained. Rather than the steep rise expected 

^ ■ within the radius of gravitational influence of a supermassive black hole, the 

^ ! random stellar velocities showed a decrease within a parsec of the center of 

K*" I the galaxy. The implied upper limit on the mass of the central black hole is 

ly-v ' only 3OOOM0, about three orders of magnitude lower than the dynamically- 

cn ' inferred mass of any other supermassive black hole. Detecting black holes of 

only a few thousand solar masses is observationally challenging but is critical 
for establishing how supermassive black holes relate to their host galaxies and 
which mechanisms influence the formation and evolution of both. 



o 



;_( ' At a distance of 850 kpc from Earth, M33 is classified (1) as a late- type ScII-III spiral, 

c/3 ■ consistent with its almost nonexistent bulge (2,3). The nucleus of M33 is very compact, 

reaching a stellar central mass density of several million solar masses per cubic parsec 

(4,5), larger than that of any globular cluster. While such high nuclear densities might 

/S ' be expected in the presence of a supermassive black hole (SMBH) (6), ground-based data 

c^ • show no evidence for a central rise in stellar velocities that would indicate the presence of 

a compact massive object in the nucleus (4). 

M33 was observed on 12 February 1999 with the Space Telescope Imaging Spectro- 
graph (STIS) on the Hubble Space Telescope (HST). Three sets of two long-slit spectra 
each, for a total exposure time of 7380 seconds, were obtained using the G750M grating 
centered on the Call absorption triplet near 8561 A (1 A corresponds to 10~^*^ meters), 
covering 19.6 km s~^ pixel"^. The pixel scale is 0''05 with a spatial resolution of 0''115 
at 8561 A. While the two spectra in each set were obtained at the same position to facili- 
tate removal of cosmic ray events, the nucleus was moved along the slit by 0''216 between 
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each consecutive set. This dithering procedure aUows for optimal correction of residual 
variations in the detector sensitivity as well as identification and removal of malfunction- 
ing pixels. The calibration steps followed the standard procedure (7) adopted for STIS 
observations of the nucleus of M32. 
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Figure 1. HST/STIS spectra of M33 extracted at the location of the nucleus 

of the galaxy (central row, marked R =0''00) and up to 0''25 (1.0 pc) off-nucleus. The 
spectra have been shifted by an arbitrary amount in the vertical direction but the same 



scale has been maintained. The decrease in the velocity dispersion of the Call absorption 
lines (8561 A) at the central positions is visible. 

The observed spectrum (Fig. 1) at every resolution element is the convolution of 
the spectra of individual stars with the line-of-sight velocity distribution (LOSVD) of the 
stellar ensemble; the latter contains information about the mean and random velocities 
of stars in the nucleus, projected along the line of sight through the galaxy. As a typical 
stellar spectrum we adopted that of the KO III giant star HD7615 which was observed 
with the same instrumental configuration. LOSVDs were extracted from the STIS spectra 
using the Maximum Penalized Likelihood (MPL) algorithm (8) and represented in terms 
of Gauss-Hermite series (9) at every slit position. The mean velocity Vq and velocity 
dispersion ao are given (modulo a standard correction applied to Vq (9)) by the first two 
terms in the Gauss-Hermite expansion. 

Neither Vq nor ao show the rise that would be expected to occur whenever a black 
hole of mass M, significantly infiuences the motion of the stars, i.e. within a distance from 
the center of the galaxy 

GM, „ f M, \ f a 

— ^ Ri 0.028" — r-^ 



r,^_^,.,^,.y__jy__^^ (1) 



as seen at the distance of M33. In particular, the central velocity dispersion is 24 ± 3 km 
s~^, significantly /ower than the dispersion of ~ 35 ± 5 km s~^ at ±0.3" ~ 1.2 pc. Given 
the ~ 0.05" resolution of STIS, equation (1) implies M, < lO^M© . 

The predicted velocities, however, depend not only on M,, but also on the gravita- 
tional potential due to the stars and on the form of the stellar orbits around the putative 
black hole (for example are the stellar orbits eccentric or circular). Therefore, a more 
rigorous upper limit to the mass of a possible black hole can be derived by constructing re- 
alistic dynamical models of the M33 nucleus, which we assume to be spherical based on its 
projected circular shape (5). While a face-on disk would also project circular isophotes, the 
position angle and ellipticity of the nucleus and the outer disk are considerably different, 
implying that the former is unlikely to be a simple extension of the latter. Furthermore, 
the properties of the M33 nucleus are consistent with those observed for globular clusters 
(4) - the prototypical spherical systems. The luminosity density (2,5,10) was represented 
by 

Kr) = ^o(-) (1 + -) . (2) 



In the above equation, a is the distance from the center of the galaxy at which the lu- 
minosity density is a fraction (in our case, one quarter) of the central value z/q, while 7 
defines the gradient in the luminosity density for r <C a. We adopt 7 = 2 and a = 1 kpc 
(5); z/Q was fixed by requiring the luminosity within 1 pc of the center to be 1.0 x lO^L©. 
The gravitational potential due to the stars was derived from z/(r) and from Poisson's 
equation as a function of the parameter M/L, the ratio of stellar mass toV— band lumi- 
nosity expressed in solar units. Past studies found M/ L < 0.5 suggestive of a young stellar 
population (4,11). 
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Figure 2. The stellar rotation curve (lower panel) and velocity dispersion profile (upper 
panel) of stars in the nucleus of M33, derived from the STIS spectra as described in the 
text. One-sigma error bars are shown. Also shown is the angular scale corresponding to 1 
parsec at the distance of M33. 



For a dynamically relaxed stellar system in a spherical gravitational potential, the 
number of stars which occupy a given location in velocity-position space, may be any 
non-negative function f{E, I) of the orbital energy E and angular momentum /. We used 
standard techniques (12) to find the / which best reproduced the kinematical data given 
M, and M/L, subject to the constraint that the integral of / over velocities reproduced 
the assumed luminosity density z/(r). The velocities predicted by / were projected onto 
the plane of the sky and convolved with the HST point spread function (PSF) and the 
STIS slit in order to allow direct comparison with the observed velocities. We emphasize 
that a two-integral (/(-E, L)) model for a spherical system, like the model adopted here, 
permits precisely as much fiexibility in the orbital distribution as a three-integral model for 
an axisymmetric system. In other words, the constraints on the mass of the M33 SMBH 
derived from a three-integral modelling code would be precisely the same as those found 
here unless the underlying stellar potential were assumed to be significantly nonspherical. 

In the absence of any additional constraints on /, we found that even black holes with 
M, > 5OOOOM0 could be made consistent with the data. However the /'s corresponding 
to these large values of M, were always found to be physically unreasonable: the stellar 
orbits changed suddenly from nearly circular at r > 0.1 pc to nearly radial at r < 0.1 pc, 
causing the projected velocity dispersion to drop sharply at a radius corresponding to the 
angular size of the STIS PSF before rising again near the black hole. After convolution 
with the PSF, the observed velocities in these solutions therefore remained low even when 
M, was large. To avoid such unphysical behavior, the solutions for / were regularized 
(12), i.e. forced to be smooth; the regularization parameter was chosen to be just large 
enough to suppress unphysical features on the scale of the PSF. 

Even after regularization, reasonable fits to the data with M, > lOOOOM© could still 
be found for values of the stellar mass to light ratio lower than 0.1 Mq/Lq. This can be 
understood qualitatively as follows. The gravitational potential is defined by the joint con- 
tribution of the stars and the central black hole. Decreasing the stellar mass-to-light ratio 
has the effect of diminishing the stellar contribution to the central potential; to compen- 
sate, the mass of the black hole needs to be increased proportionally. Decreasing M/L also 
has the more subtle effect of requiring the stellar orbits to become predominantly radial at 
large distances. As a result, the predicted line-of-sight velocity dispersion drops suddenly 
just outside of the fitted region, contrary to what is observed: the velocity dispersion in 
the M33 bulge appears to remain high, a ~ 34 km s~^, within R ^ 80" (13). By forcing 
the rms line-of-sight velocity dispersion to be greater than 30 km s~^ in the radial range 
0''5 < R < 20", such unphysical solutions were excluded. 



d 



\ o 
© 



to 
d 



Csj 

d 




logio M. (Mq) 



Figure 3. Contours of constant x^ measuring tlie fit of the dynamical models described 
in the text to the data of Figure 2. Horizontal axis is the assumed mass of a central black 
hole and vertical axis is the mass-to-light ratio of the stars in the bulge. The lowest plotted 
contour is at x^ = 1.4 and contours are separated by 0.2. 
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The best-fit / subject to these constraints was computed over a grid of {M,, M/L) 
values (Fig. 3). For M, < lOOOM© and M/L ^ 0.35, the mass of the putative black 
hole is too small to significantly affect the observable velocities and the quality of the fit 
is independent of M,. As M, is increased above ~ 2OOOM0, x^ increases as the rise in 
the central value of a provides a progressively worse fit to the data at R < 0''1. Values 
of M, as large as ~ 3OOOM0 imply x^ ^ 1.7 when only the data in the innermost two or 
three points are compared with the model, and produce a best-fit model that overpredicts 
the velocity dispersion at each of the four innermost points by ~ twice the measurement 
uncertainty in a. We take 3OOOM0 as our upper limit on M,. This value is roughly ten 
times smaller that the value ~ 5 x 10'^ Mq inferred from ground-based data (4). 

There is an empirical relation between the masses of SMBHs and the properties of 
their host bulges called the "M, — a relation," which is expressed as (14) 

1-3 f ,, / _J (3) 



IO8M0 ■ V200kms 

with a = 4.80 ±0.54. Here ac is the stellar velocity dispersion measured within an aperture 
of radius re/8 centered on the nucleus and re is the projected radius containing 1/2 of the 
light of the bulge. Because r^, is not well measured for the M33 bulge (estimates range 
from 0.5 kpc (2) to 2 kpc (10)), and the dependence of a on radius is not known accurately 
outside of the central ~ 1 pc, we conservatively adopt 21 km s~ ^ o^c ^ 34 km s~ , the 
range of values measured between 0" and 80" (13). The M,—a relation then predicts a mass 
in the range 2600 < M,/Mq < 26300, consistent with our upper limit. In contrast, the 
shallower M, —a relation proposed by Gebhardt et al. (15) would imply M, > 2.5 x 10^ M0, 
which is not consistent with the upper limit on M, obtained here. 

The M, — a relation is used to study SMBH demographics and constrain models of 
black hole formation and evolution (16,17,18). The low-mass end of the relation is of 
particular importance because SMBHs larger than ~ 10^ M0 are believed to originate 
through physical processes different from those regulating the formation of smaller mass 
black holes (19,20). However, all supermassive black holes detected so far have masses 
M, > 1O^M0 (Fig 4). Evidence for "intermediate-mass black holes" (IMBHs), with masses 
in the range 1O^M0 < M, < 1O^M0, is so far only circumstantial and relies on speculations 
concerning the nature of the super-luminous off-nuclear X-ray sources (ULXs) detected in 
a number of starburst galaxies (21,22). The connections between the possible black hole 
in M33 and ULXs is tantalizing: the M33 nucleus itself contains the brightest ULX in the 
Local Group (23), and has optical and near-infrared colors and spectra consistent with 
those of a young cluster (3) with size and mass similar to those measured for the cluster 
containing the brightest of the M82 ULXs (5,24). However, because our upper limit on 
M, in M33 is consistent with the M, — a relation as defined by much brighter galaxies, 
we cannot yet conclude that the presence of a black hole in M33 would require a different 
formation mechanism from that of the SMBHs detected in other galaxies. 
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Figure 4. The thick soHd Hne represents the M, — a relation as derived by Ferrarese 
& Merritt (14), with 1— a confidence limits on the slope shown by the dashed lines. The 
upper limit for the black hole mass in M33 (shown by the arrow) is consistent with this 
relation but inconsistent with the shallower relation advocated by Gebhardt et al. (15) 
and shown by the thin dotted line. 
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